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We examine the distribution of the magnetic anisotropy (MA) experienced by a magnetic impurity 
embedded in a metallic nano-grain. As an example of a generic magnetic impurity with partially 
filled d-shell, we study the case of d 1 impurities imbedded into ordered and disordered Au nano¬ 
grains, described in terms of a realistic band structure. Confinement of the electrons induces a 
magnetic anisotropy that is large, and can be characterized by 5 real parameters, coupling to the 
quadrupolar moments of the spin. In ordered (spherical) nano-grains, these parameters exhibit 
symmetrical structures and reflect the symmetry of the underlying lattice, while for disordered 
grains they are randomly distributed and, - for stronger disorder, - their distribution is found to 
be characterized by random matrix theory. As a result, the probability of having small magnetic 
anisotropies Kl is suppressed below a characteristic scale A e, which we predict to scale with the 
number of atoms N as Ae ~ 1/N 3 ^ 2 . This gives rise to anomalies in the specific heat and the 
susceptibility at temperatures T ~ Ae and produces distinct structures in the magnetic excitation 
spectrum of the clusters, that should be possible to detect experimentally. 


Magnetic thin films and nano-sized objects are es¬ 
sential ingredients for high-density magnetic recording. 
Magnetic nanoparticles, in particular, are considered as 
most likely building blocks for future permanent mag¬ 
nets [1-3]. Similar to molecular electronics devices [4] or 
thin metallic layers [5, 6], spin-orbit (SO) coupling plays 
an essential role in nano-particles: by restricting the free 
motion of the magnetic spins and eventually freezing 
them [7], it enables spins to store magnetic information. 
Understanding the behavior of the magnetic anisotropy 
in systems with quantum confinement is hence of crucial 
importance for nanoscale magnetic materials science. 

SO coupling induced magnetic anisotropy (MA) ap¬ 
pears to be surprisingly large in certain nanoscale and 
mesoscopic structures. A sterling example, where con¬ 
finement induced MA provides explanation for the ob¬ 
servation, is the suppression of the Kondo effect in thin 
films and wires of certain dilute magnetic alloys [5]. As 
revealed by a series of experiments on magnetically doped 
thin metallic films [5, 6], SO coupling combined with a 
geometrical confinement of the electrons’ motion induces 
a ’dead layer’ in the vicinity of the surface, where the 
motion of the otherwise free spins is blocked by MA. The 
thickness d of this ’dead layer’, consistently explained in 
terms of surface induced spin anisotropy [5, 6], depends 
on the particular host material and dopants used, but it 
can be unexpectedly large, in the range of d ~ 100A 

In confined geometries, spin-orbit (SO) interaction can 
induce magnetic anisotropies by two fundamentally dif¬ 


ferent mechanisms. In metallic compounds of heavy el¬ 
ements with strong SO interaction, the geometry of the 
sample is imprinted into the spin texture of the conduc¬ 
tion electrons’ wave function. This spin texture varies 
in space close to the surface of the sample, and induces 
a position dependent magnetic anisotropy for the mag¬ 
netic dopants. The corresponding host-induced mag¬ 
netic anisotropy (HSO mechanism), intensively studied 
in atomic-scale engineering, is presumably at work in 
magnetically doped noble metal samples, where it gives 
rise to a relatively short-ranged confinement-induced 
magnetic anisotropy close to sample surfaces [7]. Much 
stronger and longer-ranged anisotropy can, however, be 
generated by the local SO coupling at the magnetic 
dopant’s d or /-level [10, 11] in case of magnetic impuri¬ 
ties with a partially filled d or / shell, respectively [10]. 
Since, in this case, the spin of the magnetic ion is en¬ 
tangled with the orbital structure of localized / and d 
states, and couples very strongly to Friedel oscillations, 
leading to the emergence of a strong MA (LSO mech¬ 
anism). While the HSO mechanism appears to be too 
weak to explain the thin film experiments, the stronger 
and more slowly decaying anisotropy induced by the LSO 
mechanism seems to give a consistent explanation for the 
experimental observations [10-12], and appears to be the 
dominant mechanism for SO coupling induced MA in 
confined structures. 

The surface-induced MA has been thoroughly studied 
in thin films and in the vicinity of surfaces. Surpris- 
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ingly little is known, however, about the structure and 
size of confinement induced MA in nano-grains. Here 
we therefore investigate the LSO mechanism in metal¬ 
lic grains and demonstrate that symmetrically ’ordered’ 
nano-grains and nano-grains with random surfaces show 
very different behaviors. In ordered nano-grains, the MA 
constants exhibit regular structures reflecting the sym¬ 
metry of the grain. Different atomic shells of the grain 
behave very differently from the point of view of magnetic 
anisotropy, which displays Friedel-like shell-to-shell oscil¬ 
lations. Adding atoms to an ordered grain and thereby 
making its surface disordered, however, changes this pic¬ 
ture completely: in such ’disordered’ grains, the conduc¬ 
tion electron’s wave function becomes chaotic, and the 
distribution of MA parameters become gradually more 
and more random. The MA distribution is then found to 
be fairly well captured by random matrix theory, and to 
be almost independent of the magnetic ion’s position. 

In the present work, we shall demonstrate these char¬ 
acteristic properties by focusing on the simplest case of 
a magnetic impurity in a d 1 configuration embedded into 
an fee Au nano-grain host of 100-400 atoms. This model 
system captures the generic properties of most magnetic 
impurities and hosts, and allows us to study the roles of 
local and host SO interactions simultaneously. We con¬ 
struct the nano-grains by placing Au atoms on a reg¬ 
ular fee lattice starting from a central site, and then 
adding ’shells’, defined as groups of atoms that trans¬ 
form into each other under the cubic group {Oh)- We 
refer to nano-grains with only filled shells as ordered (or 
spherical) nano-grains, while nano-grains with partially 
filled outmost shells shall be referred to as disordered (or 
non-spherical) nanoparticles. We also define the core of 
the grain as the group of atoms having a complete set 
of first neighbors. To describe the electronic structure of 
Au nano-grains, we use a tight binding model with spd 
canonical orbitals and incorporate SO coupling of the 
host atoms non-perturbatively. More technical details 
can be found in the Supplemental Material [16]. 

To investigate the local SO-induced anisotropy, we 
shall use the approach of Ref. [9], and account for local 
correlations on the magnetic impurity by means of a gen¬ 
eralized Anderson model [17], which we embed into the 
Au grain described above. Similar to Anderson’s model, 
our impurity Hamiltonian (the so-called ionic model [18]) 
contains three terms: the impurity term, the conduction 
electron, and the hybridization terms (see Ref. [16]). In 
the ground state d 1 configuration, by Hund’s third rule, 
the strong local SO interaction aligns the angular mo¬ 
mentum of the d-electron antiferromagnetically with its 
spin, thus forming a D s O spin j = 3/2 multiplet. This 
multiplet remains degenerate in a perfect cubic environ¬ 
ment, and - in group theoretical terms - it transforms 
according to the four-dimensional Ts double representa¬ 
tion of the cubic point group [19]. 

Next, we need to embed this impurity into the host. 


Following Anderson, we consider hybridization of the 
deep D 3 / 2 multiplet only with s-type host electrons, 
since these latter dominate the density of states near 
the Fermi-energy (see Ref. [16]). Local cubic sym¬ 
metry implies, however, that only linear combinations 
of neighboring s orbitals, transforming as j ~ 3/2 
can hybridize with the deep Z ? 3 / 2 states. The proper 
|s 3 / 2 ), lsi/ 2 ), |s— 1 / 2 ), |s_ 3 / 2 ) basis set has been con¬ 
structed in Ref. [9, 12] and is reproduced in the Sup¬ 
plemental material, Ref. [16]. Considering then charge 
fluctuations to the d° state and performing a Coqblin- 
Schrieffer transformation [20, 21], we finally arrive at the 
following simple exchange Hamiltonian, 

'Hlso = J ^2 4n s m'\m'){m\ . ( 1 ) 

Here the {|m)| refer to the states {§,§,—§,—§} of 
the impurity, and s' m creates appropriate host electrons, 
while J denotes the strength of the effective exchange 
coupling (see also Ref. [16]). 

To handle the exchange interaction J, we can use a di¬ 
agrammatic approach similar to Ref. [9]. The dominant 
contribution to the MA is, however, simply given by the 
Hartee term, generating the effective spin Hamiltonian, 

H L = ^2 K mm'\m')(m\ 
m,m' 

with the anisotropy matrix A mm / expressed as 

/ £ F 

Pmm' ( £ ) • ( 2 ) 

-OO 

Here p L (s) denotes the local spectral function matrix of 
the symmetry adapted host operators, s' m , and Sp stands 
for the Fermi-energy. In practice, we evaluate the integral 
(2) in terms of the Green’s functions of the host [23]. 

Although disordered nano-grains do not possess spatial 
symmetries, time reversal (TR) symmetry is still present, 
and implies that, apart from an unimportant overall shift 
A'q, the anisotropy matrix K mm ' can be parametrized in 
terms of five real numbers, A M {p = 1...5), 
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which we shall refer to as the LSO-MA parameters. We 
note that the absence of SO coupling on the host atoms 
further simplifies the structure of H L , and the matrix 
elements A 3 , A 4 and A 5 vanish in case the up and down 
spin channels do not mix in the host. 

The MA matrix in Eq. (3) has two Kramers degenerate 
eigenvalues, A + = — A_, whose splitting can be used to 
define naturally the magnetic anisotropy constant as 

K L = (A+ - A_)/2 = . 


( 4 ) 
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FIG. 1: (color online) Left: Anisotropy parameters in the E-plane for an ordered grain of 225 atoms (87 core atoms). Right: 
Anisotropies in the E-plane in case of 100 disordered Au nano-grains. An ordered cluster is created by 225 atoms that fill 
completely the first six shells of the fee structure, and with 11 extra atoms, randomly placed onto the next shell. The triangular 
structure of the anisotropy parameter distribution can still be observed. 


If the magnetic impurity is placed in an ordered (spher¬ 
ical) nano-grain then the MA matrices can be different 
for different sites even if they belong to the same shell. 
Their trace, eigenvalues and, therefore, the MA constant 
should, however, be the same for all sites within the same 
shell due to the underlying (O/,) symmetry of the grain, 
as indeed confirmed by our numerical simulations, dis¬ 
cussed below (for additional information, see Ref. [16]). 

To find connections between the elements of K mm ' for 
an ordered nano-grain, we express Eq. (3) in a multipolar 
basis. Time reversal symmetry implies that K mm i can be 
expressed solely in terms of even powers of the j = 3/2 
spin operators, J x ,J y ,J z . In fact, the parameters in (3) 
couple directly to the usual 5 (normalized and traceless) 
quadrupole operators allowed by time reversal symmetry, 
Qi,...,Q 5 proportional to 2J?- Jj- J%, J%-J y , J X J Z + 
J Z J X , J x Jy + JyJ x and J y J z + J z J y , respectively. The 
local Hamiltonian can be simply expressed in terms of 
these as 

H L = ^ , (5) 

with the coefficients K t , of the Q-matrices forming a 5- 
dimensional vector. Under cubic point group transforma¬ 
tions, the first two components, ~ ( Qi,Q 2 ) and the last 
three components, ~ (Q 3 ,Q 4 ,Qs) transform into each 
other according to the E and T 2 representations of the 
cubic point group, respectively [19]. Correspondingly, 
for atoms on the same shell of an ordered grain, the 
anisotropy parameters K\ 2 transform into each other, 
and form regular patterns of triangular symmetry in 
the (Ki, A^-plane, referred to as the E-plane in what 


follows. The left panel of Fig. 1 shows the computed 
(Ki,K 2 ) values, plotted in the E-plane for an ordered 
grain of 225 atoms (87 core atoms). Throughout this 
work, we use J = 0.25 eV, a value consistent with a 
Hondo temperature below 0.1 K. Different colors denote 
the MA parameters of clusters with magnetic impurities 
placed on the different shells. Similarly, the anisotropy 
constants K 3 , K 4 . and K 5 (shown in Fig. 4 in Ref. [16]), 
induced by the SO interaction on the host Au atoms, are 
related for atoms on the same shell, and show regular pat¬ 
terns in a three dimensional space, the T 2 -space. These 
parameters are, however, smaller by about one order of 
magnitude compared to the parameters K\ 2 , implying 
that the MA constant, Eq. (4), is dominated by the E- 
type parameters. 

Next, let us examine the distribution of the magnetic 
anisotropy constants in case of disordered nanoclusters. 
First we created 100 disordered nanoclusters by adding 
11 extra atoms to an 225-atom ordered cluster, and plac¬ 
ing them randomly on the next shell of 24 possible sites. 
We then calculated the MA parameters on all core sites 
for every nanoparticle. In Fig. 1 (right) we show the 8700 
E-plane parameters obtained this way. Different colors 
represent data from different shells. Small ’clouds’ are 
observed with obvious remains of the three-fold symme¬ 
try, but there is no strictly ordered structure left anymore 
in the E-plane. 

We then increased the structural disorder of the nano¬ 
grains further, and added 25 extra atoms to an ordered 
cluster of 225 atoms, by placing them randomly on the 
next three shells. As shown in the inset of Fig. 2, for 
these strongly disordered clusters the distribution be- 
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comes almost isotropic in the -E-plane, and the trian¬ 
gular symmetry is almost entirely lost. The main panel 
of Fig. 2 shows the radial distribution of the magnetic 
anisotropy parameters, \JKf + K 2 in the E-plane. The 
observed distribution agrees very well with the predic¬ 
tions of a simple Gaussian theory, where we assume that 
the components of K E = (E'i, R' 2 ) have independent and 
Gaussian distribution, 

p(K E ) ~ e -K ® /A ® , (6) 

with the E-plane anisotropy scale A E defined as A 2 E = 
(K|;). We find that in these disordered grains the ra¬ 
dial distribution in the three-dimensional T 2 -space pa¬ 
rameters can also be fitted by a similar Gaussian en¬ 
semble, although with a smaller characteristic radius, 
VW?) =A t <A e (see Ref. [16]). The overall distribu¬ 
tion of the anisotropy K E is therefore strongly suppressed 
at small values. In the absence of host SO coupling, 
it scales as p{K E ) ~ K E for small anisotropy values, 
K l < A e , while in the presence of it p(K L ) is suppressed 
as p(Kl) ~ |El | 4 for Kl < At- This implies that typ¬ 
ical sites in a disordered grain have a finite SO-induced 
anisotropy of size ~ A E , and of random orientation, al¬ 
most independently of their precise location within the 
grain. 

We remark that, - even after adding a single extra 
atom to an ordered cluster, - the distribution of the 
eigenenergies of the host Hamiltonian agreed with the 
predictions of random matrix theory and, in agreement 
with the experimental findings [13], exhibited level re¬ 
pulsion according to a Gaussian symplectic (GS) level 
spacing distribution. The observed GS distribution re¬ 
flects the chaotic nature of the electron’s wave function 
as well as the presence of host SO coupling [16]. Building 
upon the chaotic nature of the electron’s wave function, 
one can obtain an estimate of A E , using Eq. (2) and 
assuming random plane wave conduction electron wave 
functions [22]. This yields the estimate, 

A E ~ J (Aso/ef)/E 3 / 2 , 

with N the number of lattice sites on the cluster and A so 
the SO splitting of the j = 3/2 and j = 5/2 impurity- 
levels. 

The anisotropy distribution ( 6 ) has a direct impact 
on the magnetic excitation spectrum of the nano-grains. 
Generating 50 (strongly) disordered nano-grains with 
225+40 atoms (87 core sites), we randomly chose and 
rotated 100 nano-clusters from the 4350 different sam¬ 
ples (see Ref. [16]), and determined the magnetic im¬ 
purities’ excitation spectrum averaging over the orienta¬ 
tion of them. All these 100 spectra were added together, 
and the procedure was repeated ten times. The obtained 
aggregated spectrum is shown in Fig. 3. The obtained 
spectra are grain specific (see inset), and reflect directly 



FIG. 2: (color online) Radial distribution of the magnetic 
anisotropy parameters (dots) in E-plane in case of A+=50 
samples with N=225+25 atoms. The continuous line presents 
the predictions of the Gaussian Orthogonal Ensemble. Inset: 
Distribution of (K\, K 2 ) in the E-plane for these 50 nano¬ 
grains. At this level of disorder the triangular structure is 
almost entirely lost. 


the MA energy’s distribution. The typical anisotropy val¬ 
ues, Ag « 0.57 THz, shift rapidly towards smaller values 
(GHz) with increasing grain size. 

The universal anisotropy distribution should be in¬ 
directly observable through thermodynamic quantities, 
too. In the presence of a random distribution of 
anisotropies, given by Eq. ( 6 ), we obtain a peak in the 
specific heat C(T) at T « 0.78 A E , and a low tempera¬ 
ture specific heat C(T ) ~ T 2 /A 2 E (see Ref. [16]), turning 
into a ~ T 5 anomaly for T <C At- Similarly, the coeffi¬ 
cient of the Curie susceptibility, ~ T\, should exhibit a 



FIG. 3: (color online) Orientational and disorder averaged 
excitation spectrum computed for an ensemble of 1000 ran¬ 
domly chosen and randomly oriented disordered nano-grains 
(N=225+40). The shape of the signal reflects the distribution 
of magnetic anisotropies while the small peak corresponds to 
transitions between the lowest Kramers doublets, split by the 
magnetic field. Inset: spectra of 100-100 selected grains [16]. 
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strong suppression below T ~ A E . 

Although the model discussed here has features that 
are specific, we believe that it captures many generic 
properties of magnetic impurities in a metallic grain, and 
thus allows one to draw general conclusions; For any mag¬ 
netic impurity of spin </, time reversal symmetry implies 
that the leading anisotropy term is of the form (5). In 
ordered grains, the distribution of the five parameters 
I\^ must always reflect the underlying lattice symmetry, 
and for cubic lattices, in particular, the couplings K# 
and Kj- are organized into triangular and cubic struc¬ 
tures, respectively. These parameters are expected to be¬ 
come random, and to exhibit multidimensional Gaussian 
distributions in sufficiently disordered grains. The sup¬ 
pression of the probability of having a small anisotropy, 
p(K e —» 0) = 0, as well as the predicted specific heat and 
susceptibility anomalies are also generic features, since 
they follow simply from the presence of randomly dis¬ 
tributed independent anisotropy parameters, K M . Our 
conclusions regarding the Schottky anomaly are thus gen¬ 
eral, though details of the low temperature scaling of 
C(T) may be system specific. 
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APPENDIX A - DESCRIPTION OF THE GOLD 
NANOGRAIN-HOST 

We have defined the structure of the Au nano-grains 
host of iV=100-400 atoms as follows: one can speak of 
an ordered grain when it has only filled shells around 
a center atom, while nano-particles with partially filled 
(outmost) shells have been referred to as disordered nano¬ 
grains. The shell is a group of atoms ( N s h ) on an fee lat¬ 
tice that transform into each other under the cubic group 
{Oh)- The shell-structure of a few fee nano-clusters is 
shown in Table I: an ordered grain built by, say, N =225 
atoms has 12 filled shells: a 1 center atom, 12 first-, 6 
second- ... and 24 twelfth-neighbors. The site (0,1.5,1.5) 
belongs to shell 9NNa, while the site (0.5, 0.5, 2) is on 
shell 9NNb, though they are at the same distance from 
the origin (center atom). If a disordered grain host of 
N = 236 atoms has 11 atoms in the outmost {12NN) 
shell (instead of N s h =24 atoms), then one can deal with 
C 24 configurations. In practice, we choose randomly only 
Ng grains from the big configuration space {Ns is typ¬ 
ically set as 50-100). It should be noted that one can 
put extra atoms not only into the first outmost shell but 
we never generate ’’holes” in a nano-cluster. In a given 
nano-grain the atoms those that have all the first neigh¬ 
bors are called as core atoms (denoted by N c in Table 
I). The core region is away from the surface of the nano¬ 
grain. If a nano-particle is built by N =225 atoms then 
it has N c =87 cora atoms (being in the center site and on 
six core shells). 

The electronic structure of the Au nano-grains will be 
described by a tight binding (TB) Hamiltonian. The 
model uses spd canonical orbitals, and the spin-orbit 
(SO) coupling of the host atoms is considered non- 
perturbatively. Specifically, the TB model uses (nearly) 
orthonormal basis functions which are localized at sites, 

(r I n; ao) = (r - R„ | aa) = ip a {r - R n )<j> a , (1) 

where n refers to the given site, the index a denotes the 


TABLE I: The shell structure of fee clusters. Label, a and b 
denote the shells where the atoms are at the same distance 
from the center atom but cannot be transformed into each 
other (under Oh)- N s h denotes the number of sites in a given 
shell, N is the total number of atoms and N c is the number 
of core sites in the cluster. 


Shell 

N 3h 

N 

N c 

Shell 

N sh 

N 

N c 

Center 

1 

1 

0 

16 NNa 

24 

405 

177 

INN 

12 

13 

1 

16NNb 

24 

429 

201 

2 NN 

6 

19 

1 

17 NNa 

24 

453 

225 

3 NN 

24 

43 

1 

17 NNb 

6 

459 

225 

4 NN 

12 

55 

13 

18 NN 

48 

507 

249 

5 NN 

24 

79 

19 

19 NN 

24 

531 

249 

6 NN 

8 

87 

19 

20 NN 

24 

555 

273 

7 NN 

48 

135 

43 

•21NN 

48 

603 

321 

8NN 

6 

141 

43 

22 NN 

24 

627 

321 

9 NNa 

12 

153 

55 

23 NN 

48 

675 

369 

9NNb 

24 

177 

55 

24 NN 

8 

683 

369 

10NN 

24 

201 

79 

25 NNa 

48 

731 

393 

11 NN 

24 

225 

87 

25NNb 

12 

743 

405 

12 NN 

24 

249 

87 

26 NN 

24 

767 

411 

13NNa 

48 

297 

135 

27 NN 

24 

791 

435 

13 NNb 

24 

321 

141 

28 NN 

24 

815 

459 

14NN 

48 

369 

165 

29 NN 

24 

839 

459 

15 NN 

12 

381 

177 






so-called canonical basis (real spherical harmonics), 

a = s t = 0 

Ot=Px,Py,Pz, > ( 2 ) 

Ot — dxy , d x z , hyz , d x 2_y2 , d$ z 2_ 1 t — 2 

t/v depends only on the azimuthal quantum number t 
and the spin quantum number is labeled by <r = ± g. 
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The Hamiltonian of the noble metal host is written as 


H = 


{KZo'a) = {e a 8 aa .8„.+ZH™ a , al )5 nn '+t n 0 


where the dimension of the matrix is M = 18 x N, e 0 
the so-called on-site energy parameter, 


(3) 


IS 


tjLS 
n aa,<y' cr' 


( aa\ L S | a'a') , 


(4) 


£ is the SO coupling parameter and t”’",are the hy¬ 
bridization matrix elements (or hopping integrals) be¬ 
tween the different orbitals. 

We note that on-site energies e s , s p , £d-E g and £ d-T 2g 
were used in case of all calculations. The hopping inte¬ 
grals to first- and second nearest neighbors were consid¬ 
ered. They depend only on the relative positions of the 
sites, i.e., 

C,a' = t a,a' ( R n' “ Rn) • (5) 

The numerical values for both s a and t a , a ’ can be found 
in Ref. [1]. The matrixelements of the SO coupling can 
easily be calculated with the help of following identity, 

LS=^{L + S_+L_S+)+L z S z . ( 6 ) 

The spin-orbit coupling parameter £ was determined 
from the difference of the SO-split d-resonance energies 


AEd — Ej— 5/2 — Ej = 3/2 , ( 7 ) 

as derived from self-consistent relativistic (SKKR) first- 
principles calculations [2]. This splitting is related to the 
strength of SO coupling as 

A E d ~ ^ ■ (8) 

For Au bulk we obtained £ = 0.64 eV. 

The Green’s function or resolvent operator of a nano¬ 
particle is defined as 

G(z) = G n a £ a , tJ ,(z) = (z-H)- 1 (9) 

that can be written as 


M 

<?(*) = E 


2 = 1 


K)fal 

z-Si 


( 10 ) 


where } and {|i>i)} stand for the eigenvalues and eigen- 
vectros of the Hamiltonian Eq. (3), respectively. 

Finally we define the density of states (DOS) as follows, 

n(e) = —— lim Tr (g(s + iS) — G(s — i6)\ . (11) 

2ni (5-vo V 1 

The numerically calculated values are shown in Fig. 1 
for an ordered nano-grain. The calculated Fermi-energy 
is £f = 7.4 eV, at the Fermi-energy the s contribution 
dominates the DOS. 



FIG. 1: (color online) The spd density of states (DOS) com¬ 
ponents of an ordered nano-grain (N=225) normalized to one 
atom. The calculated Fermi-energy is £f = 7.4 eV. 


APPENDIX B - DERIVATION OF THE LSO 
MODEL HAMILTONIAN 

Next we derive the local spin-orbit (LSO) model for a 
d*-type magnetic impurity. In the (non-degenerate) An¬ 
derson impurity model [3] a single energy level, e d , is con¬ 
sidered at the impurity. This level can be at most doubly 
occupied. Because of the Coulomb repulsion, the energy 
of the doubly occupied state is 2e d + U. Moreover, since 
the wavefunction of the d level is not orthogonal to the 
states of the conduction band, they may be hybridized. 
We consider a free ion that has the following four possible 
states: the d level is empty in state |d°), it is occupied 
by an electron with spin a in state Id*), and it is doubly 
occupied by electrons with opposite spins in state |d 2 ). 
The energies of these states are e(d°) = 0, e(d*) = £ d 
and e(d 2 ) = 2sd + U. The ground state of the ion has a 
magnetic moment if the magnetic doublet is the lowest 
in energy, that is if e(d 2 ) > e(d*) and e(d°) > e(d*), i.e. 
if —\U < e d + \U < 4{7. When U is large, no double 
occupancy happens, and the spin fluctuates. 

The so-called ionic model is a possible generalization 
of the Anderson model [4]. The Hamiltonian contains 
three terms as usual: the impurity , the host (conduction) 
electron, and the hybridization terms. If the multiplet of 
the impurity ion is denoted by |n, m n ), where n indicates 
the number of electrons in the shell and {m n } the set 
of quantum numbers characterizing the multiplet, and 
the corresponding energy is denoted by E n rrin , then the 
impurity term can be written as 

E E„, mn \n,m n )(n,m n \ . (12) 

n,m n 

If the impurity has only a single non-degenerate d level 
and the Coulomb interaction U is sufficiently large, so 
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that the doubly occupied configuration has high energy, 
so that the only relevant configurations are n = 0 or 
n = 1, then the Hamiltonian of the ionic model, 'Hionio 
can be approximated as 

Bo|0,0)(0,0|+X;Si, m |l,m)<l ) m|+X;^4, rn ^ )TO + 

m k,m 

+ E (El 1 ’ ™>< 0 ’ °l %m + VistjO, 0 )< 1 , m|) , 

k,m 

(13) 


find the relation between standard F 8 basis 


( l s 3/2) \ 

1 ^ 1 / 2 ) 

I 5 —1/2) 

V I s — 3 / 2 ) / 


(16) 


and the tensor-product basis {D s \a)} : where <5 = x 2 — 
y 2 or 2 z 2 — x 2 — y 2 and er =f, 4 -- We should find the 
transformation matrix between the two basis as follows, 


{\s m )} = Q{D s \a)} . (17) 


where m runs over 2 j +1 values. In case of j = \ the U = 
00 Anderson model is recovered, see Ref. [5], Sec. 1.9. 
In Eq. (13) the operator st creates a host conduction 

k,m 

electron with wavenumber k, pseudospin m and energy 
c r. Vr-s denote the s-d hybridization matrix elements. 

As we mentioned in the main paper, the desired ef¬ 
fective Hamiltonian, which describes the interaction of 
the magnetic impurity and the host (conduction) elec¬ 
trons, should be invariant under the cubic group symme¬ 
try. Therefore we have to use the F 8 symmetry adapted 
combinations of the host s orbitals (d-type or, more pre¬ 
cisely, E-type combination of these orbitals). 

The gold host atoms form an fee lattice: an impurity 
has 12 nearest neighbor host atoms (in the core region 
of a nano-grain). The operators sj, y , s' xy , s' xy , s' xy , s\. z , 

4 , 4 , 4 ) 4 z> 4 ’ 4 ) 4 create the appropriate s- 
electrons at the nearest neighbor sites with either spin 
a. The subscript of the creation operators refer to the 
nearest neighbor sites. For instance if the impurity takes 
place at site ( 0 , 0 , 0 ) then index xy denotes the coordinate 
of the site (|,|,0). Using the combinations 


s 


s 

s 


e 

xy 

e 

xz 


e 

yz 


(■4 + s xyj 


1 

E 

1 

( S i/z + 4) 


xy 

at 


\/2 


(V- + s 4 

I O x y I <J X y J 


A s xz ^ , s xz + s xz ) , 


S * z V2 


V2 

71 ( 


S- —I— o"f _ 
yz 1 yz 


(14) 


To find the proper Q, we compared the action of a rota¬ 
tion around the z and x axis by angle | in both bases 
and obtained 


Q 


/ouo 0 \ 
0 0 10 
0 0 0 -1 
V 1 0 0 0 / 


(18) 


implying 


|S3/2) = -D x 2_y2 \ |) , 
kl/ 2 ) = D 2 z2-x 2 -yA t) , 
k- 1 / 2 ) = -D 2 z'2-x 2 -y 2 \ I) 1 
k- 3 / 2 ) = D x 2 _ y 2 \ t) . (19) 


By probing for the rest of the generating point group 
elements of the cubic point-group it can easily be shown 
that this basis forms indeed a T 8 representation of the 
cubic double point-group. 

The ionic Hamiltonian, Eq. (13), without kinetic en¬ 
ergy of the host electrons reduces to 


'Hlso = E d yjm)(m\ + l^E (l TO )(°l s m + 4l°)( m l) > 

m m 

( 20 ) 

where V is the hybridization parameter, and we choose 
E o = 0 and E i >m = Ed for each m. By using a Coqblin- 
Schrieffer canonical transformation, see Refs. [ 6 ] and [5], 
Sec. 1.10, for the Hlso Hamiltonian, Eq. (20), we obtain 
that 


the d-like combinations, D s , can be expressed as 

D xz = ^ (s e xz - s% z ) , D yz = 4 ( s ^ _ s y ; (i5) 

etc. Including spin variables this leads to ten D s \a) 
combinations. The T 2 -type combinations D xz , D yz , D xy 
transform as T 5 , while the E-type combinations D x 2_ y 2 
and D 2z 2 - X 2_ y 2 transform as T 8 . As we mentioned in 
the main paper, the T 8 combinations are relevant for the 
construction of the LSO effective Hamiltonian. To con¬ 
struct the LSO symmetry adapted quantities, we should 


7-Llso = J E 


star, 


>X„ 


m,m' 


( 21 ) 


where the Hubbard operators, X m f m = |? 77 .')(m|, refer to 
the states {§, 5 ,— 5 , — §} of the impurity, and creates 
appropriate host electrons, while J denotes the exchange 
constant, 



( 22 ) 


where J was set as 0.25 eV (being consistent with the 
Hondo temperature below 0.1 K). This procedure is a 
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natural generalization of the derivation of the s-d model 
by Schrieffer and Wolff for case of j = \. 

We note that the LSO Hamiltonian Eq. (21) is ba¬ 
sically the U = oo limit of the ionic model. Consid¬ 
ering the large Coulomb interaction of the impurity d- 
level, the n = 2 double occupation is not allowed and 
the above Hamiltonian cannot be obtained in this simple 
form. However, we are convinced that a finite U can re¬ 
sult in a magnetic anisotropy (MA) of the same order of 
magnitude, i.e. the LSO model is the basic mechanism 
of the partially (but not half) filled magnetic impurities 
embedded into noble metal nanosystems. In Refs. [7] 
and [8] Au(Fe) and Cu(Fe) reduced dimensional dilute 
magnetic alloy systems (thin films) were analyzed, while 
Ref. [9] dealt with Cu(Mn) system. In case of Fe impu¬ 
rities, the suppression of the Hondo effect has been ob¬ 
served but this effect for Mn has not been measured yet. 
This observation is in agreement with our expectation 
that the LSO model can produce large enough magnetic 
anisotropy to explain the reduction of the Hondo effect, 
and suggests that the finite U does not have importance. 

APPENDIX C - CALCULATION OF THE 
MAGNETIC ANISOTROPY 

Here we derive the MA matrix from the LSO model. 
We should note that the host Hamiltonian, Eq. (3), must 
be modified in the presence of a magnetic impurity that 
breaks the two-dimensional translation symmetry. The 
simplest way to account for this constraint is to shift the 
on-site d-state energies of the impurity £ l a far below the 
valence band and add the following term to the Hamil¬ 
tonian, 

AH = A K£ a , a , = (si - s a ) 5 n0 6 n , 0 6 aa ,6 aa , , (23) 

where the impurity is at site n = 0 of a nano-grain. 

According to the LSO model, we need the Green’s 
function only for a cluster of sites, C, consisting of nearest 
neighbor atoms around the impurity and of the impurity 
itself (12+1 atoms). The corresponding Green’s function 
matrix can be evaluated as 

g(z) = g'(z) (i - A H'g(z)^ , (24) 

where I is a unit matrix and g'(z) = {G(z)}c , and G(z) 
is defined by Eq. (10). The spectral function matrix of 
cluster C is then defined as 

p c {s) = -J~- ]im (<j(e + iS) - g(s - iS)) . (25) 

27TZ (5—>0 

The dimension of the matrices defined in Eqs. (24) and 
(25) is (13x9x2), because the number of sites in the small 
cluster C is 13, and the orbital and spin-indices are 9 and 
2, respectively. Our impurity model is restricted to the 


hybridization between the impurity and the s-type con¬ 
duction electrons, hence, from the s-components of the 
pc(s) matrix we define the following projected matrix, 

p n s -c,„Ae) = Pc n s*, s Ae), (26) 

where n,n! £ C and p s -c is a (2 x 13) x (2 x 13) matrix, 
considering the the up (f) and down (J,) spin channels. 

Then we construct a 4 x 4 p*(e) matrix from p s -c re¬ 
ferring to its elements by <5i, 62 and a, a’ indices, where 
= x 2 — y 2 and 82 = 2 z 2 — x 2 — y 2 are the E-type 
orbitals. Next we transform this matrix into the r 8 sym¬ 
metry adapted basis by using Eq. (18) as follows, 

p L (e) = Qp*(e)Q' . (27) 



a - b. 

FIG. 2: First and second ordered self energy diagrams of the 
impurity spin. The dashed and continuous lines denote the 
propagators of the spin and the conduction electrons, respec¬ 
tively. 

To calculate the splitting of the four states, we employ 
Abrikosov’s pseudofermion representation [10]. Using the 
up to second order in J, the diagrams are shown in Fig. 
2, and the self-energy at T = 0 temperature is given by 

£ mro ,( w = 0) = E« ro , + E^ m ,, (28) 

where 

= J [ ^pL'(e) , (29) 

and 

/ £F POO 1 

d£ / d£ 'n— p P™m'(z)^2P™"m"(£') , 

-OO Jsf £ £ m „ 

(30) 

where (e) are the elements of p L (e) computed in the 
absence of the exchange interaction, i.e. J = 0, and £f 
is the Fermi-energy [11]. Interestingly, already the first- 
order contribution to the self-energy gives a nonvanishing 
anisotropy in the vicinity of a surface or at a site of a 
nano-grain. Therefore, as what follows we consider this 
term only. In the main paper we identify the MA matrix 
as the resulted first order self-energy, 

{K mm ,} = { E^,} , (31) 
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FIG. 3: (color online) Calculated MA constant, Kl, values for 
both ordered and disordered nano-balls. The ordered cluster 
has N=225 atoms, the solid lines correspond to the values of 
Kl calculated for the ordered case. The number of core sites 
is N c =87, i.e. the MA is for atoms located in the center site 
and on the first six core shells, see the different colors. The 
number of all sites in disordered nano-balls is A=265, i.e. 40 
extra atoms are put to the next three outmost shells. The 
number of sample is Ns = 50. 

and the effective spin Hamiltonian can be written as 

H l = Y, K mrn/ \m')(m\ . (32) 

m,m' 

The structure of H L is given by Eq. (3) in the 
main paper, the MA matrix can be parametrized by 
five real numbers, K^ (/1 = 1...5). The difference of 
its Kramers-degenerate eigenvalues defines the LSO mag¬ 
netic anisotropy constant Kl (see Eq. (4) in the main 
paper). The MA constant should be the same for sites 
in the same shell because of the Oh symmetry relations. 
The thick solid lines in Fig. 3 show the corresponding 
MA values for different (core) shells in case of an ordered 
nano-grain with 225 atoms. If one randomly puts 40 ex¬ 
tra atoms on the sites of the next three outmost shells, 
constructing disordered nano-grains, then the MA is not 
shell-degenerate anymore. 

It should be also noted that in case of tetragonal sym¬ 
metry, e.g. when the magnetic impurity is in the vicinity 
of a surface of a film or bulk material, we obtain that 
the MA has only diagonal non-zero elements, i.e. the en¬ 
ergy difference of the 3/2 and 1/2 states defines the MA 
constant (Kl = K t ). Moreover, in case of cubic sym¬ 
metry, when the magnetic impurity takes place in the 
bulk, Ki =0 and, therefore, Kl=0 as well in agreement 
with the statement that the D 3//2 ground state (T 8 ) re¬ 
mains degenerate in a cubic crystal held. This implies 
that the magnetic anisotropy should go to zero at the 
inner shells when increasing the size of nano-grains like 


in the bulk material (the characteristic energy scale of 
the MA goes by N~ 3 ^ 2 , see main paper). The calculated 
values in Fig. 3 are very fluctuating, which suggests that 
we are far from the bulk behavior. We also note that 
in case of ordered nano-clusters zero value was obtained 
for the central atom in agreement with the theoretical 
investigations based on symmetry analysis. 

The parameters of the MA matrix can be calculated for 
each core site in a given nano-grain. The matrix elements 
are different for different sites even in the same shell, but 
symmetry relations were found between the sets of the 
MA parameters. Let us enumerate the core sites by index 
i = 1...87. We label expression (32) by a site index i, i.e. 
we now introduce a set of the if/'-s. The parameter set 
{A' m } will have a site index i, too. Let i and j be two 
sites in the same shell and let the two sites be connected 
by the group element g £ Oh where gR,; = R ? . So the 
transformation rule for the matrices iT, L and can be 
written as 

H L = e i V n J H L e - lV n J ( (33) 

where J = (J x , J Vl J z ) is the angular momentum-vector 
operator, while g corresponds to the rotation around axis 
n with angle (j). From Eq. (33) we can derive relations 
between the anisotropy matrixelements of the different 
sites. 

Motivated by the quadrupole-decomposition of the 
anisotropy matrices, it seemed reasonable to calculate 
the transformation rules for the coefficients in Eq. (5) 
of the main paper. The quadrupolar operators form the 
basis of a real five-dimensional vector space (parameter 
space), and the transformation matrices between the pa¬ 
rameter sets of the different sites i and j are defined as 

K j =Y^ j \g)K\ (34) 

where is a vector formed by the set of 

The transformation matrices F^^L (gj have a very spe¬ 
cial structure, namely they are block-diagonal with two- 
and three-dimensional blocks: E^^\g) and (<?) 

in order. The five-dimensional parameter space is de¬ 
composed into two (dim=2+3) subspaces. 

Whereas the two-dimensional E parameter subspace 
has been discussed in the main paper, we focus here just 
on the three-dimensional T 2 -space. We note that the K 3 , 
K 3 and K 5 parameters span the T 2 space being by about 
one order smaller in magnitude then the parameters on 
the A-plane. Fig. 4 shows the calculated MA parameters 
in the T 2 space for the ordered nano-grain with 225 atoms 
(87 core atoms). We can identify the tetrahedral ordering 
pattern of the parameters in agreement of the structure 
of the T 2 transformation matrices. The ordered structure 
can be split even putting only one extra atom to the 
outmost shell of the spherical cluster. 

Similarly to the case of A-plane, we examine the radial 
distribution of the MA parameters in the T 2 parameter 
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FIG. 4: (color online) The magnetic anisotropy parameters 
in the T 2 -space in case of an ordered nano-grain (X=225, 
N c = 87, iVs=l). The color coding is the same as in Fig. 3. 



IK t I [meV] 

FIG. 5: (color online) Radial distribution of the magnetic 
anisotropy parameters (dots) in T 2 -space in case of TVs =50 
samples with AI=225+25 atoms ( N c =87 ). The continuous 
line presents the prediction of the Gaussian Unitary Ensem¬ 
ble. The obtained MA energy scale: At=0.13 meV. 


space. In the simulations we dealt with the same 50 
clusters as in case of Fig 3. in the main paper. The 
calculated TVspace parameters are shown in Fig. 5. The 
observed distribution agrees again with the prediction of 
Gaussian theory, similarly to Eq. (6) in the main paper. 
Vector Kt is defined as K t = (-K 3 , Ka, K-/), and the T 2 - 
space anisotropy scale At can be obtained as yj (K^) = 
At- The T 2 data in Fig. 5 can be fitted by Gaussian 
unitary ensemble (GUE), see Appendix D for a definition. 
For the (normalized) distribution we obtained 5.103 for 
the a 2 , and 1.986 for the b 2 . However we found a better 
fit in the E-plane (it was fitted by GOE), there is no 
doubt that the power is /3=2 for T 2 closer to the origin. 
At is really smaller than A e- At=0.13 meV, while At 


was 2.37 meV. 


APPENDIX D - LEVEL SPACING 
DISTRIBUTION OF THE 
HOST-EIGENENERGIES 


The spectral statistics of the Hamiltonian of the dis¬ 
ordered mesoscopic systems is usually analyzed in terms 
of random matrix theory (RMT). Here we analyze the 
energy level spacing distribution of the Hamiltonian Eq. 
(3) in the presence and in the absence of the host SO cou¬ 
pling. If £ = 0 in Eq. (3) then we speak of nonrelativistic 
grains, while by choosing the value ^ = 0.64 eV relativis¬ 
tic grains are considered. The matrix H is an M x M 
matrix, where M = 18 x N which has at most M/2 dif¬ 
ferent eigenvalues {ei} because of Kramers-degeneracy. 
The level spacings {s^} are defined as 

s i = £i +1 — Si , (35) 


and £1 < s 2 < ... < £m/ 2 - I n case of disordered nano¬ 
grains the eigenvalues are strictly doubly degenerate, so 
we always deal with M/2 different eigenvalues and, there¬ 
fore, M/2—1 level spacing values {s,}. In a given sample 
there are Ng x M/2 eigenvalues and Ng x {M/2 — 1) level 
spacings. The distribution of {e^j-s is nothing else then 
the DOS, Eq. (11), while the distribution of {s,}-s will 
be denoted by p(s). We can introduce the universal p(s) 
function as follows, 

Pp(x) = exp (- bpx 2 ) , (36) 


where the P denotes the appropriate Gaussian ensem¬ 
ble, and has a crucial role in classification of different 
universal cases, x = s/ A, where A is the expected 
value of the level spacings. In case of p = 1 one can 
speak about the Gaussian Orthogonal Ensemble (GOE), 
the P = 2 case defines the Gaussian Unitary Ensemble 
(GUE), while p = 4 corresponds to the Gaussian Sym- 
plectic Ensemble (GSE). If the time reversal symmetry 
of a system is broken, e.g. there is an external magnetic 
field, then p = 2, and the GUE can be the proper en¬ 
semble. In the presence of time reversal symmetry for 
particles with integer spins, the GOE is the appropri¬ 
ate ensemble, i.e., an orthogonal matrix diagonalizes the 
Hamiltonian. Without host SO interaction the distribu¬ 
tion of the level spacings is described by GOE. In case 
of particles with half integer spins, and in the presence 
of the host SO coupling, one has to deal with GSE [12]. 
We note that the constraints 

/»oo /»oo 

/ pp(x)dx = 1 and / xpp(x)dx = 1 (37) 

J 0 J 0 

imply that 


re+\(p + 2)/2) 

r*+i(08 + i)/2) 


and bp 


F 2 ((/3 + 2)/2) 
F 2 ((/3 + l)/2) ’ 
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where F is the gamma function, r(n) = (n — 1 )!, and 


7r 7r 

(ii = 1.571, h = - ; 

32 4 

02 = — ~ 3.242, 62 = — 

71 ' 7T 


_ 262144 

4 7297T 3 


11.597, 64 


0.785 (GOE) , 
i 1.273 (GUE) , 

= ^ m 2.264 (GSE) . (39) 


The calculated density of states, n(e), in case of a nano¬ 
grain host of iV=225 Au atoms can be seen in Fig. 1. 
The large DOS arises from the ’d-band’ around 3 eV. 
Let us denote the integrated DOS by N(e). Calculating 
N(e) we observed two energy ranges, £d and £ Sl where 
the integrated DOS is a nearly linear function and it can 
be supposed, that n(e) is constant in these ranges, at 
least on averages. The average level spacing can then be 
defined as, 


A.,; = 


1 




(i = d , s). 


(40) 


In case of d-band between 3 and 4.5 eV ( £d ) a good sta¬ 
tistical analysis can be done because of the large values 
of n(£d)- The other energy range, £ s is chosen from 6 to 
14 eV, where the dominant component of the DOS is the 
s-type states. The value of n(£ s ) is smaller than n(£d ), 
therefore we cannot do a satisfactory statistical analysis. 
Importantly, the range £ s contains the Fermi energy. 

Kubo-relation says that the expected value of the level 
spacings is inverse proportional to L d , where L stands for 
the linear size of a d dimensional system. It follows that 

A N,i ~ ; (i = d, s) . (41) 

Reassuringly, the ratios of, say, A 393 ^ = 2.337 meV and 
A 273 y d = 3.667 meV agree with the value of ||| within an 
accuracy of 10 % in case of nonrelativistic particles and 
the ratio of A 237 y d = 3.983 meV es Ai 6 5 ; d = 5.469 meV 
agree with within an accuracy of 5 % in the presence 
of host SO coupling. 

The distribution of the level spacings is shown in Fig. 
6 . We dealt with 225+25 atoms in a given cluster and 
chose randomly 25 sites from the next 96 sites at the 
12NN, 12NNa and 13NNb shells (see Table I), which 
creates Cff possible configurations. In our simulation 
Ns = 50 relativistic samples were calculated. The data 
are well fitted by GSE: a4=19.13, 64 = 2 .867 with 95 % 
confidence bounds. 

It should be noted that the largest deviation of the 
calculated data from the universal p(x) function can be 
observed in the tail region. The deviation of the host level 
spacing statistics (fat tail) from the prediction of RMT 
together with the anomalies observed in anisotropy pa¬ 
rameter space (angle correlations) strongly suggest that 
the effect of the cubic symmetry of the underlying lattice 
is not negligible compared with total randomness. Most 



x, normalized level-spacing 


FIG. 6: (color online) The level spacing distribution for 
As=50 disordered nano-grains with A=225+25 atoms in the 
d-band. The average level spacing is A = 3.9 meV. The data 
are well fitted by GSE function, 04=19.13, 64 = 2 .867 with 95 
% confidence bound, 0.965, and are in agreement with the 
experiments [13]. 


importantly, as we mentioned in the main paper, the ob¬ 
tained level spacing distribution in Fig. 6 is in agreement 
with the experiments [13]. 


APPENDIX E - MAGNETIC EXCITATION 
SPECTRUM OF IMPURITIES IN NANO-GRAINS 

Here we calculate the excitation spectra of magnetic 
impurities embedded in random nano-grains. The LSO 
MA energy can be calculated due to the effective Hamil¬ 
tonian H l , see Eqs. (2) and (5) in the main paper. To 
describe the spectra of an arbitrarily oriented sample we 
rotate H L by an angle ip about the unit vector n. Then 
the Hamiltonian becomes 

H^ v = e ^* H L e -iv nJ (42) 

where J = (J x , J y , J z ) is the angular momentum-vector 
operator. We solved the eigenproblem of this Hamilto¬ 
nian H^ v \i) = Ei\i) (i = 1,2,3,4). The Fermi’s golden 
rule for the intensities is given as 

S{u) oc ^2 S(uj - (.Ej - Ei)/h)\(j\J x \i)\ 2 pi, (43) 

j>i 

~P E i 

where pi = g— is the Boltzmann weight of the *th 
eigenstate, and Z = e~ l3Ei is the partition function 
with /3 = 1 /kgT the inverse temperature. 

We simulated the excitation spectra of disordered 
nano-grains with A=225+40 atoms according to the fol¬ 
lowing procedure: we generated A(g=50 nano-grains, and 
for each such grain we placed one magnetic ion on the dif¬ 
ferent core sites. There are 87 such sites, which means we 
generated 4350 different samples. From these 4350 nano¬ 
grains we have randomly chosen 100 clusters, and rotated 










them with angles <p uniformly distributed in the interval 
[0,7r] and chose the vector n uniformly distributed on the 
unit sphere. The Dirac-delta in Eq. (43) was replaced by 
an (unnormalized) Lorentzian with half-width T = 0.1 
THz for the absorption spectra (for better comparison 
with experiments we give the frequencies in THz). The 
temperature was set to T = 4.2 K. All these 100 spec¬ 
tra were added up, and the procedure was repeated ten 
times. The obtained results are shown in Fig. 4 of the 
main paper. 


APPENDIX F - SCHOTTKY ANOMALOUS 
HEAT CAPACITY FOR NANO-GRAINS WITH 
MAGNETIC IMPURITIES 

Here we derive the specific heat (heat capacity) of the 
nanoballs with magnetic impurities. We have seen that 
the MA matrix has two Kramers degenerate eigenstates 
with the energy splitting equal to 2AA, where Kl is the 
MA constant, see Eqs. (3-4) in the main paper. The 
canonical partition function for a single impurity can be 
written as 


Z = 2(1 + e- 2KhP ) , (44) 


where f3=(ksT) 1 is the inverse temperature. The ther¬ 
modynamic value of the energy is defined as 


E(T) 


dlnZ _ 2 K l 

d/3 1 + e 2KL P 


(45) 


Fig. 2 in the main paper demonstrates that the main (E- 
plane) contribution to the MA constant follows a GOE- 
type random distribution, 

p(K L ) = aK L e~ bK i . (46) 


with Ae = see Eq. (6) in the main paper. Note 
that a=2b because p(e) is normalized to unity. Hence 
Eq. (46) can be written as 


P(K L ) = -£2 K lZ ■ (47) 

E 

By using Eq. (45), we calculate the average of the energy 
over this distribution, 

{E(T)) = jf E^+-+K L )iK L . (48) 

Inserting Eq. (46) into Eq. (48), and introducing the new 
variable x = Kl/3, we obtain for the expectation value of 
the energy, 


(E(T)} 


4k%T 3 

^E 



2 

are 


0 2x 


-dx 


(49) 


or, alternatively, 


/■oo 2 -a 2 x 2 

<- E(a )) = 4A B a 3 / — 2x dx , (50) 

Jo 1 + e 

where we introduced the normalized temperature a = 
pjppr- The averaged specific heat, 


(C(T)) 


d(E(T)) 
dT 


can then be expressed as a function of a as 


(51) 


<C(a)> 



8ck 2 (E 4 ) 


1 + e 


2x 


dx . 


(52) 


To get the temperature dependent specific heat, the in¬ 
tegration in Eq. (52) has to be performed numerically 
for different values of a. The result of this procedure can 
be seen in Fig. 7. Evidently, the random (GOE) distri¬ 
bution of the MA constants induces a Shottky peak in 
the specific heat at a « 0.78, i.e. at T m » 0.78 A e- It 
can also be noted that C(a ) is proportional to a 2 , thus 
C(T) ~ T 2 / A\ for temperatures, T <C T m . 



FIG. 7: (color online) Specific heat (in units of fcs) of nano¬ 
grains with magnetic impurities, obtained from the numerical 
integration in Eq. (52) as a function of the normalized tem¬ 
perature, a = ksT/AE- This curve is universal for all dis¬ 
ordered nano-grains, where the MA constants follow a GOE- 
type of distribution, Eq. (47). The grain-specific informations 
are hidden in the parameter A®. 
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